Equations of motion, initial and boundary conditions for GRB 
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We compare and contrast the different approaciies to tiie optically thick adiabatic phase of GRB 
all the way to the transparency. Special attention is given to the role of the rate equation to be self 
consistently solved with the relativistic hydrodynamic equations. The works of Shemi and Piran 
Piraii, Shemi and Narayan Meszaros, Laguna and Rees [3| and Ruffini, Salmonson, Wilson and 
Xue are compared and contrasted. The role of the baryonic loading in these three treatments 

is pointed out. Constraints on initial conditions for the fireball produced by electro-magnetic black 
hole are obtained. 
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I. INTRODUCTION 

The phenomenon of gamma ray bursts (GRBs) at- 
tracts a lot of attention in the modern relativistic as- 
trophysics. In spite of revolutionary progress in obser- 
vational technics and instruments over the past decades, 
GRBs remain misterious events since their discovery. Al- 
though the nature of the energy which powers GRBs re- 
mains in general uncertain, there exists a conventional 
model which describes the phenomenon since the mo- 
ment when the huge energy is released in a compact spa- 
tial region. It is referred to as the cosmic fireball model 

a. 

The qualitative study of cosmic fireballs was given in 
the pioneering paper by Cavallo and Rees Here we 
remind briefly their main results. The authors considered 
sudden energy release in a compact spherical region of 
space in the form of photons. If the initial average photon 
energy in units of electron rest mass energy in the fireball 
is larger than unity. 
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the fireball is opaque due to large number of electron- 
positron pairs produced. These pairs in turn produce 
new photons, sharing their energy, so the systems cools 
at the expense of growing number of particles. It be- 
comes transparent only when e < 1, and pairs annihi- 
late. However, this behavior changes if admixture of a 
plasma is also present in the fireball. This admixture is 
conveniently parametrized by 
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where ctt is the Thomson cross section, Rq is initial 
radius of the fireball and n , n'^as initial number den- 
sities of electon-positron pairs and ionized gas (which we 
will call simply plasma in the following) respectively. 

The plasma is assumed to be in the form of ionized 
hydrogen and electons, totally neutral, and it could be 
enough plasma to make the fireball optically thick even 
after the energy decreases below the pair production 
threshold. Then the fireball starts to expand adiabat- 
ically and its energy diminishes as inverse radius until 
the condition t ~ 1 is fulfilled. Once the optical depth 
reaches unity, the fireball becomes transparent and pho- 
tons stream away. This is a basic qualitative model, 
which was developed in great details in subsequent years. 
Unfortunately it was abondoned later in favor of other 
models implementing for instance the notion of beaming. 
In particular, great attention was paid to the interaction 
of expanding relativistic baryons left from the fireball 
with interstellar medium (ISM) surrounding the fireball. 

The aim of the paper is to give a critical review of ex- 
isting models for isotropic relativistic fireballs, compare 
and contrast these models. The paper is organized as 
follows. In the next section we give the basic equations 
and describe approximations involved into description of 
the fireballs. Then we present our model which differs 
from models in the literature as it describes the dynamics 
of the fireball taking into account the rate equations for 
electron-positron pairs. Then we compare various mod- 
els for the fireball. In the final section we discuss initial 
conditions for the fireball and impose corresponding con- 
straints arising from our model. Conclusions follow. 



where Eq,M,c are initial total energy of the fireball, 
total mass of the plasma and the speed of light. So initial 
optical depth of the fireball is given by 



II. ENERGY-MOMENTUM PRINCIPLE 

The basis of description for relativistic fireballs is the 
energy-momentum principle. It allows to obtain rela- 
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tivistic hydrodynamical equations, or equations of mo- 
tion for the fireball, energy and momentum conservation 
equations which are used extensively to describe inter- 
action of relativistic baryons of the fireball with the in- 
terstellar gas, and boundary conditions which are used 
to understand shock waves propagation in the decelerat- 
ing baryons and in the outer medium. Consider energy- 
momentum conservation in the most general form[19i]: 



(4) 



where Fj^^ are Cristoffel symbols and g is determinant 
of the metric tensor. Integrating over the whole three- 
dimensional volume we obtain 



T^"'.^dV = 0. 



(5) 



Integrating over the whole four-dimensional volume and 
applying divergence theorem we get |8j] 

T^^'.^dVdt = / T'"'\^dV = 0, (6) 
tJv ' Jv 

where are covariant components of the outward drawn 
normal to the three-dimensional hypersurface (volume 
V). 

Define the momentum four- vector P^: 
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From ([7]) and ([5]) in Minkowski metric (when F(^^ = 0) 
we see that 



dP° 

dt 
dP^ 

~dr 



dt 



dV = - 



dV = - 



(10) 



so if the energy and momentum fluxes through the sur- 
face S bounding considered volume V are absent the en- 
ergy and momentum are constants during system evolu- 
tion. Supposing this is the case we arrive to the conser- 
vation of energy and momentum: 



P^ = const. 



(11) 



This equation is important to describe interaction of the 
baryons left from the fireball with the interstellar gas. 
Assume the energy-momentum tensor in the form of the 
ideal fluid 



(12) 



where lo = e + p is proper entalpy, p is proper pressure 
and e is proper energy densities. Now suppose spherical 



symmetry [20], which is usually done for fireballs descrip- 
tion. Using spherical coordinates with the interval 



ds^ = -df + dr^ + r'^-de'^ 
we rewrite ([4]): 
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dt 



j,2 Qj. 



{r^T^^) - r (P33 + p44 sin2 6) = 0, (15) 



arriving to equations of motion for relativistic fireballs 

dij'^uj) dp 1 d , 2 2 \ n 

+ I2fll(^ H =0' (16) 



dt dt dr 
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dp 



dt 
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dr 
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where the four- velocity and the relativistic gamma factor 
are defined as follows: 

C/^ = (7,7^^,0,0), 7=(l-7.2)-i/2, (18) 

the radial velocity u is measured in units of speed of light 
u — v/c. 

Now suppose that there is a discontinuity on the fluid 
flow. Suppose the three-dimensional volume is a spher- 
ical shell and choose the coordinate system where the 
discontinuity is at rest so that in ^ for normal vectors 
to the discontinuity hypersurface Aq. we have 



A„A" = 1, 



Ao = 0. 



(19) 



Let the radius of the shell Rg be very large and shell 
thickness A be very small. With — > oo and A — » 
from (O we arrive to 



(20) 



where the brackets mean that the quantity inside is the 
same on both sides of the discontinuity surface. This 
equation together with continuity conditon for particle 
density flux [nW] — was used by Taub Q to obtain rel- 
ativistic Rankine-Hugoniot equations, describing shock 
waves dynamics. Relativistic shocks are supposed to ap- 
pear during collision of the baryonic material left from 
the firebaU with the ISM [g. The origin of the afterglow 
could be connected to the conversion of kinetic energy 
into radiative energy in these shocks [1, [ll|, [l^, 13 ■ How- 
ever, our scenario differs from that, namely we suppose 
that fully radiative condition during this interaction is 
satisfied. Our model allows to explain the afterglow phe- 
nomenon without consideration of shocks [14] as sources 
of radiation. 



III. QUASI-ANALYTIC MODEL OF GRBS 

The first detailed models for relativistic fireballs were 
suggested in the beginning of nineties [ll_3i 9- Inde- 
pendent calculations performed in Q and|5| give precise 
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understanding and we describe our approach first, men- 
tioning the main differences with the existing literature. 

First of all, the source of energy, being obscure in previ- 
ous models, is supposed to be the energy extraction pro- 
cess from electromagnetic black hole (EMBH) J^] . The 
second difference is that initially not photons but pairs 
are created by overcritical electric field of the EMBH and 
these pairs produce photons later. This plasma, referred 
to as pair-electro-magnetic (PEM) pulse expands initially 
into the vacuum surrounding the black hole reaching very 
soon relativistic velocities. Then collision with the bary- 
onic remnant of the EMBH takes place and the PEM 
pulse becomes pair-electro-magnetic-baryonic (PEMB) 
pulse (see lj| for details). It was shown that the fi- 



nal gamma factor does not depend on the distance to 
the baryonic remnant and parameters of the black hole. 
The only crucial parameters are again the energy of the 
dyadospherc, or simply Eq, and baryonic admixture: 



B = 



= V 



(21) 



The exact model is based on numerical integration 
of relativistic energy-momentum conservation equations 
(|16ll7p together with the baryonic number conservation 
equation [21 



(22) 



However, the most important distinct point from all 
previous models is that the rate equation for electron- 
positron pairs is added to the model and integrated si- 
multaneously in order to reach self-consistency. 

Here we concentrate on the simple quasi-analytical 
treatment presented in 0,11] (see also [lj|). The PEMB 
pulse is supposed to contain finite number of shells each 
with flat density profile. The dynamics is determined by 
the following set of equations: 
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(23) 
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where F ~ | is a thermal index giving the pressure p — 
(F — l)e, e is a proper internal energy density e = p — pB, 
Pb — nBTUpC^ is baryon mass density in comoving frame, 
V and V are the proper volume in the comoving frame 
and the volume in the coordinate frame: V — j{r)V. 
Subscript "0" denotes initial values, and all quantities are 
assumed being averaged over finite distibution of shells 
with constant width and density profiles. All components 
such as photons, electrons, positrons and plasma ions give 
contribution to energy density and pressure. This set of 



equations is equivalent to ([32]) and ((34)) (see below). The 
next step is to take into account the rate equation: 



dt 



-N,. 



), (26) 



where a is the mean pair annihilation-creation cross 
section, is the thermal velocity of e^-pairs. The co- 
ordinate number density of e^-pairs in equilibrium is 
Ng± (T) — 7ng± (T) and the coordinate number density 



of e*-pairs is N^± 



7ne±. 



For T > irieC^ we have 



n^±{T) = n-y{T), i.e. the number densities of pairs and 
photons are equal. The pair number densities are given 
by appropriate Fermi integrals with zero chemical poten- 
tial, at the equilibrium temperature T. 

For an infinitesimal expansion of the coordinate vol- 
ume from Vo to V in the coordinate time interval t — to 
one can discretize the last differential equation for nu- 
merical computations. 

The most importants outcomes from analysis per- 
formed in are the following: 

- the appropriate model for geometry of expanding 
fireball (PEM-pulse) is given by the constant width 
approximation (this conclusion is achieved by com- 
paring results obtained using (|16ll7p and simplified 
treatment described above), 

- there is a bound on parameter B which comes from 
violation of constant width approximation, B < 
10-2 (?7 > 10^). 

The last conclusion is crucial since it shows that there 
is a critical loading of baryons, when their presence pro- 
duce a turbulence in the outflow from the fireball, its 
motion becomes very complicated and the fireball evolu- 
tion does not lead in general to the GRB. 

Exactly because of this reason, the optically thick fire- 
ball never reaches such large radius as = fo?/^ (dis- 
cussed in 0, see section (IVII) ) since to do this, the bary- 
onic fraction should overcome the critical value Be — 
10^2. For larger values of Be the theory reviewed here 
does not apply. This means in particular, that all con- 
clusions in [3| obtained for r > Vb are invalid. In fact, for 
B < Be the gamma factor even does not reach saturation 
(see sec. IVip . 

The fundamental results coming from this model are 
the diagrams presented at fig. [T] and [51 The first one 
shows basically which portion of initial energy is emitted 
in the form of gamma rays Ej when the fireball reaches 
transparency condition r ~ 1 and how much energy gets 
converted into the kinetic form of the baryons Ek left 
after pairs annihilation and photons escape. 

The second one gives the value of gamma factor at the 
moment when the system reaches transparency. 

The energy conservation holds, namely 
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(27) 
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FIG. 1: Relative energy release in the form of photons emitted 
at transparency point E^/Eq (solid line) and kinetic energy 
of the plasma E^/Eq (dashed line) of the baryons in terms 
of initial energy of the fireball depending on parameter B — 
rj~^ obtained on the basis of quasi-analytic model. Thick 
line denotes the total energy of the system in terms of initial 
energy Eq. 



45 mp 1 1 
0.032, 



(29) 
(30) 



where nip, me, are proton and electron masses, go — 
a — dimensionless temperature T and radius TZ 

2 t 

of the fireball are measured in units of and Ae = — ^ 
correspondingly, and the subscript " 0" denotes initial val- 
ues. The temperature at transparency point in the case 
when plasma admixture is unimportant is nearly a con- 
stant for a range of parameters of interest and it nearly 
equals 



Tp = 15 keV. (31) 
Adiabatic expansion of the fireball implies: 
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where £ — 



is a radiative energy in terms of electron 



rest-mass energy. From the energy conservation ([7]) , sup- 
posing the fluid to be pressureless and its energy density 
profile is constant we have in the coordinate frame: 
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FIG. 2: Relativistic gamma factor of the fireball when it 
reaches trasparency depending on the value of parameter B. 
Dashed line gives asymptotic value 7 — B~^ . 



const. 



(33) 



Supposing at initial moment 70 = 1 and remembering 
that Etot = E + Mc^ we arrive to the following funda- 
mental expression of relativistic gamma factor 7 at trans- 
parency point: 



Clearly, if the baryons abundance is low most energy is 
emitted when the fireball gets transparent. It is remark- 
able that almost all initial energy is converted into kinetic 
energy of baryons already in the region of validity of con- 
stant thickness approximation B < 10^^, so the region 
10~^ < B < 10^^ is the most interesting from this point 
of view. 



£0 + Mc^ 



£ + Mc^ (%f)^ + l 



(34) 



M 



where A4 

One can use this relation to get such important char- 
acteristics of the GRB as observed temperature and ob- 
served energy. In fact, they can be expressed as follows: 



IV. SHEMI AND PIRAN MODEL 

In this section we discuss the model, proposed by 
Shemi and Piran This quantitative model gives 

rather good general picture of relativistic fireballs. 

Shemi and Piran found that the temperature at which 
the fireball becomes optically thin is determined as 

Tesc=mmiTg,Tp), (28) 

where Tg and Tp is the temperature when it reaches 
transparency with respect to gas (plasma) or pairs: 



%bs = "fTesc, (35) 
T b 

£obs = ^OTf — • (36) 
-'esc 

These results are presented at fig. |3l In the limit of 
small T] we have 7 = (1 -I- 77), while, for very large t] 
the value of gamma factor at transparency point is 7 = 
To/Teac, and it has a maximum at intermediate values of 
77. We donote by dashed thick line the limiting value of 
7] parameter rjc = B~^ . For ij < rjc the approximations 
used to construct the model do not hold. It is clear that 
because of the presence of bound rjc the value "f = r] can 
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FIG. 3: The relativistic gamma factor (upper dashed line), the 
observed temperature (solid line), and the ratio of observed 
energy to the initial energy of the fireball (lower dashed line) 
as a function of rj (see [3| ) ■ The values of parameters are the 
same as in the cited paper. Thick dashed line denotes the 
limiting value r^c- The values of rj when gamma factor reaches 
maximum and gets constant are also shown. 



be reached only as asymptotic one. In effect, the value 
r]c cuts the region where saturation of the gamma factor 
happens before the moment when the fireball becomes 
transparent. 

It was found that for relatively large ij > 10^ the pho- 
tons emitted when the fireball becomes transparent carry 
most of the initial energy. However, since the observed 
temperature in GRBs is smaller than initial temperature 
of the fireball, one may suppose that a large part of initial 
energy is converted to kinetic energy of the plasma. 



V. SHEMI, PIRAN AND NARAYAN MODEL 

Piran, Shemi and Narayan Q present generalization of 
this model to arbitrary initial density profile of the fire- 
ball. These authors performed numerical integrations of 
coupled energy-momentum relativistic consevation equa- 
tions (|16ll7p and baryon number conservation equation 
(I22p . They were mainly interested in the dependence of 
the observed temperature, gamma factor and other quan- 
tities on the radius of the fireball. Their study results in 
the number of important conclusions, namely: 

- the expanding fireball has two basic phases: a ra- 
diation dominated phase and a matter-dominated 
phase. In the former, the gamma factor grows lin- 
early as the radius of the fireball: 7 oc r, while 
in the latter the gamma factor reaches asymptotic 
value 7 ~ ?7 + 1 . 

- the numerical solutions are reproduced with a good 
accuracy by frozen-pulse approximation, when the 
pulse width is given by initial radius of the fireball. 



FIG. 4: The ratios of radiation and kinetic energy to the 
initial energy of the fireball predicted by Meszaros, Laguna 
and Rees model. Thick line denotes the total energy of the 
system in terms of initial energy. Energy conservation does 
not hold. 




FIG. 5: Relative energy release in the form of photons emitted 
at transparency point Eots/Eo of the GRB in terms of initial 
energy of the fireball depending on parameter B — rj~^ . Thick 
line represents numerical results and it is the same as in fig. 
[T] Normal line shows results for the analytic model of Shemi 
and Piran ^ . Dashed line shows the difference between exact 
numerical and approximate analytical results. 



(37) 



where A ~ i?o is the width of the leading shell with con- 
sant energy density profile, R is the radius of the fireball. 

They also present the following scaling solution: 



R = 
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(38) 
(39) 



The last conclusion is important, since the volume V 
of the fireball can be calculated as 



where subscript "0" denotes some initial time when 
> few, which can be inverted to give "f{R)- 
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VI. MESZAROS, LAGUNA AND REES MODEL 

The next step in developing this model was made in 
Q. In order to reconcile the model with observations, 
these authors proposed a generalization to anisotropic 
(jet) case. Nevertheless, their analytic results apply to 
the case of homogeneous isotropic fireballs and we will 
follow their analytical isotropic model in this section. 

Starting from the same point as Shemi and Piran, con- 
sider ((32)) and pi]) . The part of the paper, containing 
analytical results, describes the geometry of the fireball, 
gamma factor behavior and the final energy balance be- 
tween radiation and kinetic energy. Magnetic field effects 
are also considered, but we are not interested in this part 
here. 

Three basic regimes are found in Q for evolution of the 
fireball. In two first regimes there is a correspondence 
between the analysis in the paper and results of 0, so 
the constant thickness approximation holds. It is claimed 
in 0, that when the radius of the fireball reaches very 
large values such as Rb = RqT]^ the noticable departure 
from constant width of the fireball occurs. However, it is 
important to note, that the fireball becomes transparent 
much earlier and this effect never becomes important (see 
section nil)) . 

The crucial quantity, presented in the paper is - 
the maximum possible bulk Lorentz factor achievable for 
a given initial radiation energy Eq deposited within a 
given initial radius Rq: 
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(40) 
(41) 



where Sq is initial baryon (plasma) mass surface den- 
sity. 

All subsequent calculations in the paper [3| involves 
this quantity. It is evident from (HDl) that the linear de- 
pendence between the gamma factor F and parameter 77 
is assumed. However, this is certainly not true as can be 
seen from fig. [31 We will come back to this point in the 
following section. 

Another important quantity is given in this paper, 
namely 



This means, above rjp the presence of baryons in the 
fireball is insufficient to keep the fireball opaque aft^er 
pairs are annihilated and almost all initial energy de- 
posited in the fireball is emitted immediately. The esti- 
mate of the final radiation to kinetic energy ratio, made 
in 3] is incorrect, because kinetic and radiation energies 
do not sum up to initial energy of the fireball thus vio- 
lating energy conservation (P7)) . This is illustrated at fig. 
|4l The correct analytic diagram is presented in fig. [5] 
instead. 



VII. APPROXIMATE RESULTS 

All models for isotropic fireballs are based on the fol- 
lowing points: 

1. Flat space-time, 

2. Relativistic energy- momentum principle, 

3. baryonic number conservation. 

Anthough the EMBH model starts with Reissner- 
Nordstrom geometry, the numerical code is written for 
the case of fiat space-time simply because curved space- 
time effects becomes insignificant soon after the fireball 
reachs relativistic expansion velocities. The presence of 
rate equation in the model has a deep physical 

ground and its luck in the other treatments means incom- 
pleteness of these models. Indeed, the number density 
of pairs influences the speed of expansion of the fireball. 
However, in this section we neglect the rate equation and 
discuss the common points between all considered mod- 
els. 

First of all, let us come back to fig. [31 For almost all 
values of parameter 77 the gamma factor is determined 
by gas (i.e. plasma or baryons) admixture according to 
(Uni), consider this case below. For given initial energy 
and radius this temperature depends only on r] only, so 
one can write: 



1] + 1 
arp + i 



(44) 



where 



r 

'-P ~ rp ■ 
J-P 



(42) 



a = 2.1 • 10^V^7^o 



0.5 



(45) 



This is just the asymptotic behavior of the gamma fac- 
tor at fig. [3] for very large 77. Using it, the authors cal- 
culate the value of 77 parameter above which the pairs 
dominated regime occurs: 



p3 

m 

p 



(43) 



From this formula we can get immediately the two 
asymptotic regimes, namely: 
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77+1, 77 < 77max, 



(46) 



Notice, that the constant a is extremely small number, 
so that after obtaining precise value of ?7max by equating 



7 



to zero the derivative of function (|44p one can expand in 
Taylor series the result and get in the lowest order in a, 
that: 



(47) 
(48) 



2\3 

a 



In particular, in the case shown in fig. [3] one has 
Vmax = 2.8 • 10^, 7rnax = 9.3 • 10** whilc according to 
(SOI) Tm = Vmax ~ 1-75 • 10^. Clearly, our result is much 
more accurate. Actually, the value r,„ in (|40l) is ob- 
tained from equating asymptotes in (|46|) and there exists 
the following relation: 
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FIG. 6: Relativistic gamma factor when transparency is 
reached. The thick line denotes exact numerical results, the 
normal line corresponds to analytical estimate from Shemi 
and Piran model, the dashed line denotes the asymptotic 
value 7 = . 



a = (rory)-!/^. (49) 

Now we can explain why the observed temperature 
(and consequently the observed energy) does not depend 
on 1] in the region rjmax < V < Vp- From the second line 
in it follows that the gamma factor in this region 
behaves as 7 ex r]~^^^, while Tesc oc 77^/^. These two ex- 
actly compensate each other leading to independence of 
the observed quantities on rj in this region. This remains 
the same for rj > rjp also, since here T^sc = Tp =const 
and from (|44| 7 =const. 

VIII. SIGNIFICANCE OF THE RATE 
EQUATION 

The rate equation describes the number densities evo- 
lution for electrons and positrons. In analytic models 
it is supposed that pairs are annihilated instantly when 
transparency condition is fulfilled. Moreover, the dynam- 
ics of expansion is influenced by the electron-positron en- 
ergy density as can be seen from Therefore, it 
is important to make clear whether neglect of the rate 
equation is a crude approximation or not. 

Using eq. p4|) one can obtain analytic dependence of 
the energy emitted at transparency point on parameter 
B and we compare it at fig. [5l 

We also show the difference between numerical results 
based on integration of eqs. ((23)) - ((26)) and analytic results 
from Shemi and Piran model. The values of parameters 
are: fj, = 10^ and ^ = 0.1 (which correspond to Eq = 2.87- 
10^'* ergs and Rq = 1.08 • 10^ cm). One can see that the 
difference peaks at intermediate values of B. The crucial 
deviations however appear for large B, where analytical 
predictions for observed energy are about two orders of 
magnitude smaller than the numerical ones. This is due 
to the difference in predictions of the radius of the fireball 
at transparency moment. In fact, the analytical model 
overestimates this value at about two orders of magnitude 
for B = 10~^. So for large B with correct treatment 



of pairs dynamics the fireball gets transparent at earlier 
moments comparing to the analytical treatment. 

At the same time, the difference between numerical 
and analytical results for gamma factor is significant for 
small B as illustrated at fig. [HI While both results co- 
incide for B > 10""^ there is a constant difference for 
the range of values 10~^ < B < 10^^ and asymptotic 
constant values for the gamma factor are also different. 
Besides, this asymptotic behavior takes place for larger 
values of B in disagreement with analytical expectations. 
Thus the acceleration of the fireball for small B is larger 
if one accounts for pairs dynamics. 

It is clear that the error coming from neglect of the 
rate equation is significant. This implies that simple an- 
alytic model of Shemi and Piran gives only qualititive 
picture of the fireball evolution and in order to get cor- 
rect description of the fireball one cannot neglect the rate 
equation. 

Moreover, the difference between exact numerical 
model QilHl ^'^d approximate analytical models [l| be- 
comes apparent in various physical aspects, namely in 
predictions of the radius of the shell when it reaches 
transparency, the gamma factor at transparency and the 
ratio between the energy released in the form of photons 
and the one converted into kinetic form. The last point 
is crucial. It is assumed in the literature, that the whole 
initial energy of the fireball gets converted into kinetic 
energy of the shell during adiabatic expansion. Indeed, 
taking typical value of parameter B as 10^"^ we find that 
according to Shemi and Piran model we have only 0.2% of 
initial energy left in the form of photons. However, exact 
numerical computations give 3.7% for the energy of 

photons radiated when the fireball reaches transparency, 
which is a significant value and it cannot be neglected. 

To summurize the above discussion, we present the re- 
sult of this survey in the TablelH It is important to notice 
again that comparing to simplified analytic treatment, 
accounting for the rate of change of electron-positron 
pairs densities gives quantitatively different results on 
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RufRni et. al. 


Shemi, Piran 


Piran, Shemi, Narayan 


Meszaros, Laguna, Rees 


conservation: 










energy- momentum , 


yes 


yes 


yes 


yes 


baryon number, 


yes 


do not consider 


yes 


yes 


rate equation 


yes 


no 


no 


no 


constant width 


justify 


do not consider 


justify 


in part 


approximation 










model for 7(r) 


numerical/analytic 


no 


numerical / analytic 


numerical/analytic 


model for 7(77) 


numerical 


analytic 


do not consider 


analytic 



TABLE I: Comparison of different models for fireballs. 



the ratio of kinetic versus photon energies produced in 
the GRB and the gamma factor at transparency mo- 
ment, which in turn leads to different afterglow prop- 
erties. Therefore, although analytical models presented 
in sections |TV] and IVIII agree and give correct qualita- 
tive description of the fireball, one should use numerical 
approach described in sec. IIIII in order to compare the 
theory and observations. 



100000. 



1000 
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IX. CAVALLO AND REES DIAGRAM AND 
EMBH MODEL 

As already mentioned in the section IIIII the fireball 
can be produced by EMBH if the electric field strength 
exceeds some critical value Ec- Within the shell sur- 
rounding just formed EMBH where E{r) > Ec called 
dyadosphere [l^, [l^ the overcritical electric field pro- 
duces electron-positron pairs which in turn produce pho- 
tons by reactions 7 + 7 <-> + e~ , quickly thermalize 
and begin to expand being accelerated by the radiative 
pressure leading finally to observed GRB. 

The fireball can clearly emerge in the EMBH theory 
only if the dyadosphere radius r^ya is larger than the 
horizon radius r+ . These quantities are defined as follows 



and 



^dya 



Qn 



(51) 



where G is Newton's gravitational constant, /i and ^ 
are parameters of EMBH, fi = M/Mq, ^ — Q/Qmax, 
Qmax = M^/G, Mq is solar mass and the critical electric 
field strength £c is defined [13, Ull as 



he 



(52) 



where e is electron charge and h is Planck's constant. 
Condition rdya = i"+ gives [15] : 



/z = 5.31 10^ 



(i + yr^)' 



(53) 



0.1 



0.01 



1 



1. xlO 0.0001 
FIG. 7: Bound on parameter space within EMBH model. 



It defines the maximal allowed mass for given charge of 
EMBH for which the electric field outside the horizon is 
overcritical one and can produce electron-positron pairs. 
The lower bound on EMBH mass comes from the maxi- 
mal mass of the neutron star M = 3.2Mq. Equality ((53|) 
can be inverted to get the lower bound on EMBH charge 
for this value, so finally we have the following ranges of 
basic parameters: 



3.2 — l-I'min ^ < /i 
2.41 10-^ = Crmn < ^ < 6na 



„,a. - 5.31 10^ (54) 

1. (55) 



(50) These bounds together with ((53)l arc illustrated at fig. [T] 

Cavallo and Rees suggested to describe the fireball 
with two parameters 9* and e, where 



Eq (Tt 



(56) 



where nip is proton mass and e is defined in ([T]). Pa- 
rameter 9* is tightly connected with optical depths with 
respect to electron-positron pairs Tpp and electrons from 
the plasma admixture Tgas- 



(57) 

nip 

9* = mgas. (58) 

The diagram presented at fig. 1 of allows to predict 
the evolution of the fireball. It is of great interest to find 



9 




1 10- 10^ 10^ 10* 



FIG. 8: Cavallo and Rees diagram for EMBH model (details 
see in the text). 



which restrictions on this picture are predicted by EMBH 
model. These constraints are shown in fig. [5] 

The vertical dashed line denotes the pairs production 
threshold so to the right from this line pairs production 
occurs within the fireball. The inclined dashed line de- 
notes the condition Tpp — 1 and the fireball with pa- 
rameters in the region above this curve is opaque due 
to pairs. The solid horizontal line denotes the condition 
Tgas = 1 and the fireball is opaque due to plasma admix- 
ture above this line. The region of parameters allowed by 
EMBH model is represented by grey region. The upper 
boundary is given by the condition ^ = Cmax while the 
lower boundary corresponds to ^ — fimin- The maximal 
possible value of parameter 6* — 3.39 10^® is shown by 
the dotted line. The maximal possible average particle 
energy produced by EMBH is e = 306. At the same time, 
admissible initial conditions with small 9* (up to 6* = 0) 
lie in a very thin region close to the pairs production 
threshold e = 1. 

The values of parameter 6* for EMBH case are much 
higher than ones considered by Cavallo and Rees. It im- 
plies that the energy is released in a much compact region 
of space and with much larger intensity. Clearly the fire- 
ball produced by EMBH cannot have e < 1, otherwise 
pairs would not be created. On this diagram the solid 
horizontal line corresponds to B = 10~^, the critical case 
for EMBH theory. For smaller values of B this line moves 



up because it represents the condition 9* = rj = B~^ . It 
can even cross the allowed region of parameters for the 
fireball for very small B. This means that below this line 
the fireball will be opaque only to pairs and not to the 
plasma, since the plasma admixture is tiny. 

From fig. [H] it can be seen that the fireball produced 
by EMBH is almost always opaque both to pairs and to 
plasma, i.e. its initial conditions lay in the the region 
III defined in [7] . This implies that the fireball necessary 
undergoes adiabatic expansion with 9* (x e^. Therefore, 
cases I and IV described by Cavallo and Rees are irrel- 
evant for EMBH model. The case II is only partially 
relevant for almost pure fireballs when the plasma ad- 
mixture is very small B <C 10~^°. 



X. CONCLUSIONS 

We compared existing isotropic models of GRBs, so 
called fireball models. It is shown that the crucial differ- 
ence between our approach and other models in the liter- 
ature is the presence of the rate equation which accounts 
for electron-positron pairs densities evolution during ex- 
pansion of the fireball. This results in quantitative dif- 
ference between predictions of our quasi-analytic model 
and analytic models in the literature. Considering its 
significance we conclude that in order to compare the- 
ory and observations it is necessary to take into account 
rate equation together with energy and mass conserva- 
tion conditions. 

Another important difference is the presence of bound 
on baryonic loading parameter Be = 10~^ which comes 
from violation of constant thickness approximation used 
in our quasi-analytic model. The same bound should 
be present in all analytic models in the literature. As a 
consequence, the broadening of the relativistic shell re- 
sulting from the fireball never happens before it reaches 
transparency. Besides, the gamma factor does not reach 
saturation and the value j = i] = B^^ is only asymptotic 
one for B < 10~^. 

Bounds on dimensionless quantities describing initial 
conditions for the fireball are deduced from EMBH the- 
ory. It is shown that the fireball appearing as the result 
of energy extraction process from electro-magnetic black 
hole is always opaque both to its electron-positron pairs 
and electons of the baryonic material which the fireball 
absorbs during early stages of expansion. As a conse- 
quence the fireball resulting from EMBH necessarly un- 
dergoes adiabatic expansion all the way to transparency. 
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